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In this letter we consider the dynamic behaviors of spin-orbit coupled Bose condensates realized 
in recent experiments. We show that there exists an interaction induced ac Hall response which 
is absent in a non-interacting system. This condensate has two distinct equilibrium phases known 
as the plane wave phase and the stripe phase. In the plane wave phase, we show that an ac 
longitudinal current will induce an ac radial current in the transverse direction, and vice versa, as 
a cooperation effect of spin- velocity locking and spin-dependent interaction. In the stripe phase, we 
show that the dominant longitudinal response to a transverse radial current is sliding of the density 
stripe, because it is the low-lying excitation mode originated from spontaneous spatial translational 
symmetry breaking in this phase. 



The Hall effect is known as the production of current 
transverse to the voltage difference, or vice versa. Con- 
ventionally, it is caused by the Lorenz force for particles 
in a magnetic field or the Coriolis force in a rotating 
frame, which couples the particle's motion in one direc- 
tion to its motion in the transverse direction. In systems 
of neutral atomic gas with synthetic magnetic field pQ [5] , 
a superfluid Hall effect can manifest itself in the collec- 
tive oscillations and has also been observed in a recent 
experiment [3]- 

With similar Raman coupling technique, an artificial 
spin-orbit (SO) coupling can also been generated in ul- 
tracold bosons [IHS] and fermions [7]. As an analogy to 
condensed matter systems, conventional Hall and spin 
Hall effects have also been predicated for some SO cou- 
pled cold atom systems [5HTT]. However, since the config- 
uration of SO coupling realized in current experiments is 
a special one as an equal weight combination of Rashba- 
type and Dresselhaus-type. Only in one spatial direction, 
where two Raman beams are counter-propagating, the 
motions of atoms are coupled to their spins. Thus, its sin- 
gle particle Hamiltonian along three different directions 
are separable and commute with each other. That means 
the motions along different directions are completely un- 
correlated if no interactions. Therefore the conventional 
Hall and spin Hall responses are absent here. 

Nevertheless, realizing SO coupling in cold atom sys- 
tems brings several new ingredients that are absent in 
condensed matter systems. First, the system can be 
bosonic; secondly, the interactions between atoms are 
spin-dependent; and thirdly, there exists a harmonic con- 
finement potential. These ingredients do give rise to in- 
triguing physics for equilibrium physics, such as stripe 
superfluid phase [TH [13] and half vortex phase P^HTT] . 
In this letter we focus on whether they will also give rise 
to non-trival dynamics which has not been discovered 
in condensed matter systems before, and indeed we find 
that there exists an alternative and more striking Hall 
response arisen from spin-dependent interactions. 

The geometry under consideration is also slightly dif- 
ferent from conventional ones. Conventional geometry of 



Hall response is shown in Fig. [ija), while in this system 
the Hall response is that an ac longitudinal current will 
induce an ac radial current in the transverse plane, as 
shown in Fig. [ijb). For a trapped cloud, this is mani- 
fested as the coupling between the radial breathing mode 
and the longitudinal dynamics. 

This dynamic behavior also strongly depends on the 
underlying equilibrium quantum phases and reveals their 
fundamental properties. It is known that there are two 
distinct superfluid phases in this system depending on ex- 
ternal parameters, which are known the plane wave phase 
and the stripe phase [HITnillH], as shown in Fig. [l|c). In 
the plane wave phase, we show that an ac longitudinal 
current couples to transverse motion. This comes from 
that the SO coupling locks spin to the velocity. While 
in the stripe phase, the dominant response is the cou- 
pling between the transverse current and the longitudi- 
nal sliding mode of density stripe. This reveals the spon- 
taneous spatial translational symmetry breaking in the 
stripe phase and the sliding mode is a low-lying Gold- 
stone mode in this phase. 

Model. We illustrate the physics using F — 1 con- 
densate as an example. Three spin states of = 1 are 
coupled by two counter-propagating laser beams along x 
direction [IHB] , and the Hamiltonian of this system is 
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V^t(r)HoV(r)+5on'(r)+g252(r) dh (1) 



where V'^(r) = (^Li(r), Vj(r), (r)), Ho = + 
'Ho,± + Vtrap, and no,x is given by 
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T~l-o,± ~ tfkj_/{2m)I, where _L denotes {y,z) plane and 
I is 3 X 3 identity matrix. Here 2hkQ is the momen- 
tum transfer during the Raman process, S is two-photon 
Raman detuning, e is the quadratic Zeeman energy. 
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FIG. 1: (a) Configuration of conventional Hall effect; (b) Con- 
figuration of Hall response proposed in this work; (c) Phase 
diagram for the system considered in this work. In PI and 
P2 regimes the single particle dispersion has two local min- 
ima, while in P3 regime there is only one minimum. The grey 
regime is the stripe phase while its outside regime is the plane 
wave phase. Q, and 5 are in unit of 2_Er — h^kQ/m. In all 
plots of this work, we take 172 = —go/W- 

n(r) = tp'^{r)ip{r) and S(r) — ■(/'^(r)S?/;(r), where S is 
3x3 Pauli matrix vector. The energy dispersion for 
the lowest eigen-state of Ho,x is denoted by e(fcj.) and 
is shown in the inset of Fig. [ijc). e{kj.) has two local 
minima for small S and small f2, and displays only one 
minimum for either large S or large fl. We consider a har- 
monic trap of cylindrical symmetry with ujy — — uJi_ 
and Vtiap = {muJxx'^ + muj±r'^) /2. Thus, the single parti- 
cle Hamiltonian is separable in three different directions 
and they commute with each other. 

Variational Wave Function Approach. A variational 
wave function approach has been successfully used for 
studying dynamics of a BEC in absence of SO coupling 
[T^ . Here we generalize it to the case with SO coupling. 



The variational wave function is taken as a product of an 
overall Gaussian profile and the spin wave function (j). 

Tj—x^y^x ^ 'I ^ 

where A = {Ax, Ay, Az) is the centcr-of-mass position, 
^tUt^) denote the width of condensate, and 5^ 
is introduced as conjugate variable to i?^ In the 

regime where e{kx) has two minima at k±, (j) is assumed 
to be 

(/) = cos6le'¥e'P+("-A)^(p+ ^)+sin6'e-*¥e»p-('--A)^(^_ ^) 

(4) 

where p± are momenta around energy minima (fc±, 0, 0), 
respectively. While in the regime where e{kx) has a single 
minimum, saying at is assumed to be 

^^g,p+(r-A)^(p^^^) (5) 

£.{p±,x) is a three-component normalized spin wave func- 
tion taken as the lowest eigenstate of Hamiltonian Eq. 

This is the main assumption we make here, that 
is, the spin wave function is always locked to the center- 
of-mass velocity during the whole dynamic process due 
to SO coupling. This assumption is supported by recent 
experimental observations 6j. 

To compute the energy with wave function Eq. Q , we 
note typically {p+^x ~ P-,x)Rx ^ Ij thus the terms that 
contain e~*^^'+ ="~''-'^)^-'^^ will be ignored. This means that 
to a very good approximation the energy functional main- 
tains spatial translational invariance, because the poten- 
tial gradient due to the harmonic trap is much weaker 
compared with SO coupling. Under this approximation, 
the condensate energy £ is given by 
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where f{d,p+^x,P-,x) denotes the interaction energy den- 
sity of a uniform system 

f{0,P+,x,P-,x)= yd3r(go|0|" + 52|rS0p) (7) 

For the equilibrium state, obviously one has Sr/ ~ 0, 
A,! = for ?7 = x, y, z and p±,j_ ~ 0. p±,x, Q and i?,, 
should be determined by coupled equations d£ /dp± — 0, 
d£/dRn = and d£/dO = 0. If the solution gives 
< < 7r/2, the condensate wave function is a super- 
position which gives rise to the stripe order, and this 
phase is named as "stripe superfluid" phase [T^l US] , as 



indicated by the grey area of Fig. [ijc). Note that / is in- 
dependent of ip due to the translational symmetry, which 
means that ip of an equilibrium state can take any value 
determined by spontaneous symmetry breaking. If the 
solution gives 9 — oy 6 = 7r/2, or if e{kx) has a single 
minimum as wave function Eq. ([5]), the system will be in 
a "plane wave superfluid" phase. The main part of this 
paper below is to use this variational approach to study 
collective modes in these two phases, respectively. 

Plane Wave Phase: For the plane wave phase, without 
loss of generality, we can take 6* = in Eq. Q in the two 
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FIG. 2: Hall response in the plane wave phase, (a) Frequency 
of three breathing modes (solid lines labelled by (a), (b) and 
(c)) and the frequency of dipole mode (dots) as a function 
of ijJx/ou± in absence of spin-dependent interaction g2. The 
arrow indicates the position of mode resonance, (b): Hall 
effect nearby resonance point, (bl): A longitudinal current 
(black solid lines) induced by excitation of transverse radial 
current (red dashed line); (b2): A transverse radial current 
(red dashed line) induced by excitation of longitudinal cur- 
rent (black solid line). J is in unit of Jo = hko/m. t is in 
unit of to = h/2Eji, and = h^kl/{2m). For this plot 
hujx = 0.64_Er, and fko = Q.2Eyk.. For instance, in recent *^Rb 
experiment, E-g.. = 2 x lO'^Hz and to is 0.04ms. 

minima regime or equivalently take Eq. ([s]) in the single 
minimum regime, and the energy £ is given by Eq. ([6]) 
with 9 ~ Q, and automatically / depends on p+.x only. 
Considering the Lagrangian C = J (Pv^* {id / dt — H)^, 
we arrive at 
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(8) 



The equation-of-motion for dipole motion along x direc- 
tion is given by 



A, 



de{p 



+ ,x) 



df/dp+.x 



-ujIAx. (9) 



And the radial dynamics of breathing mode is given by 



Rrq — Rrj5y^ 



TT^R^RxRyRz 



(10) 



The / function contains both go term and g2 term. 
In the plane wave phase, the density is a constant, and 



thus the go term is also a constant. The momentum de- 
pendence of / comes from the 52 term, because the spin 
configuration changes with momentum due to SO cou- 
pling and thus the spin-dependent interaction energy also 
changes. If (72 is much smaller than 170 ^nd one can ignore 
52 term, then df /dp+^x is zero, the dipole mode Eq. ^ 
and the breathing mode Eq. ( 10 1 are completely decou- 



pled. In this case, the breathing mode is the same as that 
in a scalar condensation derived in Ref . [12] , and the fre- 
quencies of three different breathing modes are shown in 
Fig. [2ja) using solid lines. While in contrast to a con- 
ventional condensate, the Kohn theorem will be violated 
due to SO coupling and the dipole mode frequency will 
be different from Wx as shown in Fig. [2][a) using solid 
lines with dots [B] . In this case, the coupled equation Eq. 
([9]) with df /dp+^x — predicts how the dipole mode fre- 
quency depends on S and il, which is very consistent with 
recent experiment on ^^Rb [6], where (?2 is only 0.46% of 
50- 

The coupling between dipole and breathing modes will 
become significant with larger 172, for instance, if similar 
experiment is done with magnetic atoms like Cr |20| , Dy 
[21 1 and Er [52]. Nevertheless, even for smaller (72, since 
there exists a window of aspect ratio ujx/^± in which the 
frequency of bare dipole mode is at resonance with one 
of the breathing modes, as indicated by arrow in Fig. 
[2ja), significant coupling between these two modes can 
be obtained near the resonance regime. In this case, the 
full response behavior requires numerically solving the 
coupled differential equations Eq. (|9| and Eq. (10). The 
results are illustrated in Fig. [2|b). 

Fig. illustrates the coupling between Jx — Ax 

and Jj^ = =^ r f$*^$-$^$*V In Fig. bl we 

2ni J \ dr± dr± J ^ 

show when we excite a radial ac current J± by suddenly 
compressing the condensate in the transverse direction, 
Jl will induce a longitudinal current Jx- Vice versa, 
in Fig. b2, we show when we excite a longitudinal ac 
current Jx by suddenly displacing the condensate away 
from the center of harmonic trap, Jx will induce a radial 
current J±. This is precisely the Hall response induced 
by interactions. We also notice that there are two dif- 
ferent frequencies in the dynamics. This is because in 
resonant regime, two nearly degenerate modes are split 
by spin-dependent interaction 92 term. Thus the dy- 
namics exhibits similar behavior of quantum beats. The 
larger frequency (smaller period) is nearly bare frequency, 
while the smaller frequency (larger period) is mostly de- 
termined by coupling provided by 92- 

Stripe Phase: Now we turn to the dynamics in the 
stripe phase. Here the dynamic equations involve both 
A± and p± . The Lagrangian is given by 



R^S,, 



008 26*1 -f. 



(11) 



where Pc.r^ = cos^ (^P+,q + sin^ Op^ ^i- Thus, the equation 
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FIG. 3: Hall response in the stripe phase, (a) The relative 
phase evolves as a function of t after a sudden radial com- 
pression in the transverse direction. The inset illustrates the 
longitudinal spatial density profile Px = J d^r^|$p at differ- 
ent time Ti, T2 and T3, which demonstrates that the time 
evolution of if corresponds to sliding of density stripe, x 
is in unit of 1/fco. (b) Time dynamics of radial current J± 
in the transverse direction after a sudden change of 0. t is 
in unit of fi/2£'R, and Eb. = fi^fco/(2m). In both case, we 
choose parameters 57 = CSOBr, 5 — — 0.24i5R, hwx = 0.33i?R 
and fkj± = 0.2iJR which gives 6 — 0.85, koR^ = 12.6 and 
k(iR±_ — 20.7 for equilibrium phase. For (a), we suddenly 
compress R± to its half value, and for (b), we suddenly change 
8 to 1.05. 



for dipole mode is modified as 
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dp^ 



The equation for breathing mode is the same as Eq. ( 10 1. 
Moreover, for the stripe phase since < < 7r/2 and 
sin 20 7^ 0, there is one more dynamic equation for the 
relative phase ip given by 



if : 



1 d£ 
sin 26 m 



' P—.x)Ax 



(13) 



Because the energy is independent of tp, there is no dy- 
namics for 6. 

Here we discuss two types of Hall response shown in 
Fig. [3] Note that for an equilibrium state, 9 is given 
by the balance between the single particle energy and 
the interaction energy. First, in Fig. [3^, if one ex- 
cites a breathing mode along the transverse direction. 



this changes the interaction energy, and therefore d£ /d9 
also becomes non-zero. The direct consequence is that ip 
will change as a function of time according to Eq. ( [l3| . 
This leads to the linear dependence of (p on time t, as 
shown in Fig. [3^. Moreover, the breathing mode also 
weakly couples to center-of-mass motion in the longitu- 



dinal direction A^, thus, the second term in Eq. 13 leads 
to small wiggles on top of the linear behavior. However, 
this effect is much weaker. Time evolution of ip means 
that stripe slides along the longitudinal direction, as one 
can see from the inset of Fig. |3^ for density profile at dif- 
ferent times. Since ip corresponds to the phonon degree 
of density stripe, which is the most low-lying excitation, 
and therefore it is the most dominative response in the 
stripe phase. Secondly, in Fig. [3]d, if we suddenly change 
9, it strongly alters the interaction energy and thus ex- 
cites radial current in the transverse direction. 

Summary: In summary, we have proposed two new 
types of Hall response in SO coupled condensates recently 
realized experimentally. These effects are remarkable be- 
cause they are purely caused by interactions: (i) for the 
plane wave phase, because SO coupling locks spin to the 
velocity, the spin-dependent interaction energy changes 
as longitudinal velocity varies, and the change of interac- 
tion energy excites breathing motion along the transverse 
direction; and (ii) for the stripe phase, because the stripe 
order comes from the superposition of two wave functions 
with different momenta, and the weight of superposition 
sin0(or cos 9) is determined by the interaction energy, 
thus the change of interaction energy by compressing the 
transverse confinement excites the dynamics of its con- 
jugate variable (p, which describes the phonon mode of 
the stripe. As one can see, the first effect is directly 
connected to spin-velocity locking and the second one is 
connected to the nature of stripe phase, thus, although 
we use current experimental setup as an example to il- 
lustrate these two effects, they may also exist in a wide 
range of condensed matter and cold atom systems with 
SO coupling or stripe order. 
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